Role of electronic excitations in magneto- Raman 
spectra of graphene 



Oleksiy Kashuba 

Institute for Theory of Statistical Physics, RWTH Aachen, 52056 Aachen, Germany 
E-mail: kashubaOphysik . rwth-aachen . de 

Vladimir I. Fal'ko 

Department of Physics, Lancaster University, Lancaster, LAI 4YB, UK 

Abstract. We investigate the signature of the low-energy electronic excitations in 
the Raman spectrum of monolayer and bilayer graphenes. The dominant contribution 
to the Raman spectra is due to the interband electron-hole pairs, which belong to 
the irreducible representation A 2 of the point group Cq v of the graphene lattice, and 
are characterised by crossed polarisation of incoming and outgoing photons. At high 
magnetic fields, this is manifested by the excitation of electron-hole (e-h) inter-Landau- 
level transitions with selection rule n~ — > n + . Weaker Raman-active inter-Landau- level 
modes also exist. One of those has a selection rule similar to the infrared absorption 
process, n~ — > (n ± 1)+, but the created e-h excitation belongs to the irreducible 
representation E2 (rather than Ei) and couples to the optical phonon mode, thus 
undergoing an anticrossing with the optical phonon G-line in Raman in strong magnetic 
field. The fine structure of acquired by the G-line due to such anticrossing depends on 
the carrier density, inhomogeneity of doping, and presence of inhomogeneous strain in 
the sample. 



PACS numbers: 73.63.Bd, 71.70.Di, 73.43. Cd, 81.05.Uw 
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1. Introduction 

Experimental observation of electronic excitations in graphene is mostly associated with 
angle-resolved photoemission spectroscopy [U [2] and optical absorption in monolayers [31 
SI El El El El M and bilayers P El HOI HH H21 H31 HH H51 US] where Raman spectroscopy 
has been used to study graphene phonons p2], and became the method of choice for 
determining the number of atomic layers in graphitic flakes [TTJ, [HI [191 EQl [211 1221 [23] . 
In particular, single- and multiple-phonon-emission lines in the Raman spectrum of 
graphene and the influence of the electron-phonon coupling on the phonon spectrum 
have been investigated in great detail [2H 123 123 [2HI I2S1 M ED 

This manuscript presents a theory of inelastic light scattering in the visible range 
of photon energies accompanied by electronic excitations in mono- and bilayer graphene 
both with and without an external magnetic field. We classify the relevant modes 
according to their symmetry, analyse peculiar selection rules for the Raman-active 
excitations of electrons between Landau levels in graphene at strong magnetic fields, 
review the fine structure of the phonon-related Raman under conditions of magneto- 
phonon resonance (when the inter-Landau- level transition is in resonance with phonon), 
and model the dependence of the spectra on the value and homogeneity of carrier density 
and strain in the sample. 

Graphene is a gapless semiconductor [33l EH EHl [36] , with an almost linear Dirac- 
type spectrum e = avp in monolayer case and parabolic e = ap 2 /2m for bilayer, in the 
conduction (a = +) and valence (a = — ) band, which touch each other in the corners of 
the hexagonal Brillouin zone, usually called valleys. In bilayer graphene there are two 
additional parabolic band with the gap 271. The conduct ion- valence band degeneracy 
in the Brillouin zone corners of the band structure of graphene is prescribed by the 
hexagonal symmetry Cq v of its honeycomb lattice, and it is natural to relate Raman- 
active modes to the irreducible representations of the point group Cq v . 

We argue that the dominant electronic modes generated by inelastic scattering of 
photons with energy Q less than the bandwidth of graphene are superpositions of the 
interband electron-hole pairs which have symmetry of the representation A2 of the group 
Cq v and are odd with respect to the inversion of time. Their excitation process may 
go in two ways both of which consist of the same steps, but in reversed order. In the 
first process the absorption of a photon with energy Q transfers an electron from an 
occupied state in the valence band into a virtual state in the conduction band, followed 
by emission of the second photon with energy Q = Q — u. In the second process the 
reversed order of events leads the the virtual electron state with large deficit of energy 
as opposed to the first process with large excess. Net amplitude is determined by the 
sum of partial amplitudes of the opposite processes. The dominance of such process over 
the one involving the contact interaction [371 EH EH] of an electron with two photons is 
a peculiarity of the Dirac-type electrons in graphene. 

Strong magnetic field leads to the quantisation of the electronic spectrum into 
Landau levels [40J. In graphene Landau levels (LL) appear symmetrically in the 
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Table 1. Classification of the Hamiltonian components by Cq v irreducible 
representations. 

conduction bands e" (a = ±) with peculiar zero-energy LLs in monolayer (n = 0) 
and bilayer (n = 0.1). For the inelastic light scattering in graphene, we find peculiar 
selection rules, n~ — > n + of the dominant Raman-active transitions (solid line in Fig. [3]), 
in contrast to the An = ±1 transitions between Landau levels which are dominant in 
the absorption of left and right-handed circularly polarised infrared photons [7J. The 
weaker Raman mode with a similar inter-LL selection rule is also differs from absorption 
process, having different symmetry properties which allow the interaction with Raman 
scattering of light on optical phonons creating anticrossings at the G-peak line, see 
Fig. |4j Raman spectroscopy, therefore, provides data supplementary to that obtained 
in optical absorption. 

The following theory is based upon the effective low-energy Hamiltonian derived 
from the tight-binding model of electron states in graphene. It describes electrons in 
the conduction and valence bands around the Brillouin zone (BZ) corners K (£ = +) 
and K' (£ = — ) [41, 42J. For monolayer graphene the effective Hamiltonian reads, 

H m = &<r ■ p - ^(o x (p 2 x - p 2 y ) - 2a v p xPy ). (1) 
070 

Here er = (o x , o y ) are Pauli matrices acting on the sublattice components of electronic 
states on A and B sublattices of graphene honeycomb lattice, and p being the in- 
plane momentum counted from the BZ corner. The first term in Eq. ([T]) determines 
Dirac spectrum electrons. The monolayer hopping parameter 70 ~ 3eV determines 
the bandwidth, ~ 670, and Dirac velocity of electrons v ~ 10 8 cm/s [TQl [TJl [191 S3]- 
The second term takes into account weak trigonal warping [41] breaking continuous 
rotational symmetry of Dirac Hamiltonian. The basis is constructed of wavefunctions 
corresponding to atomic sites A, B in the valley K, and B, A in K' . Since 4-spinors 
{(pAKi <Pbk, <Pbk', PAK'} realise 4D irreducible representation of the full symmetry 
group of the crystal, matrix operators can be combined into irreducible representations 
HU [261 [27] of the group C 6v , see Table (T) 

Bilayer graphene, Fig. [T]d, consists of two coupled sheets of monolayer with AB 
(Bernal) stacking (as in bulk graphite [45]), and has symmetry D3/, which differs from 
Cq v by accompaniment of a z — > —z reflection to the rotations by tt/3, comparing to 
monolayer graphene, but, otherwise, D 3 h and C% v are isomorphic. Its unit cell contains 
four inequivalent atoms Al, Bl, A2 and B2 where numbers denote the layer. The 
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Figure 1. Schematic of the monolayer (a) and bilayer (b) graphene crystal lattice, (c) 
The band structure of mono- and bilayer graphene in the vicinity of the K point along 
the p x axis. Also shown are some of the two-step processes leading to the creation 
of an electron-hole pair in the low-energy bands accompanied by the absorption of a 
photon followed by emission. Note that the energy of any intermediate state (white 
circle) is much smaller than photons energy 17, Vt. 

conventional Hamiltonian in the vicinity of the valleys is 

* - ( " ^ \ ) , (2) 

where, t> 3 ~ O.lf is related to the weak direct Al — > B2 interlayer hops j3U [TjJJ 03], 
and 7i ~ 0.4eV is the direct interlayer coupling [H El El [10l [121 UHl SH SSI S3]. The 
basis is constructed using components corresponding to atomic sites Al, B2, A2, Bl in 
the valley K and B2, Al, Bl, A2 in K'. 

The electronic dispersions around the K point for mono- and bilayer are compared 
in Fig. [I] Two low-energy bands touch each other at the centre of the valley, which is 
also the position of the Fermi energy in the neutral structure (a neutrality point), and 
are described by the effective two-band low-energy Hamiltonian written in the basis of 
orbitals on the sites Al and B2 [36] . 

„,2 



H, 



beS 



7i 



y(pl-P 2 y ) + 2o-yp xPy ). 



(3) 



Two other, referred to as high-energy bands, are split by the interlayer coupling 71 from 
the neutrality point. 



2. Inelastic light scattering in graphene at zero magnetic field 

To include the interaction of the electrons with photons in the effective Hamiltonians ((!]) 
and pi), we construct the canonical momentum p — - A, where the vector potential of 
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Figure 2. Feynman diagrams for Raman scattering amplitude R from ground state 
to final state with electron-hole excitation. 
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includes an annihilation operator &q l3z ,i for a photon with in-plane momentum q, energy 
(which determines its out-of-plane momentum component q z = \JVl 2 /c 2 — q 2 ), and 
polarisation 1. Expanding the Hamiltonians up to the second order in the vector 
potential one obtains the interaction part, 

o2 



H: 



int 



c 



e" 
2t^ 
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(e/c)d Pi H is the current vertex and ^sdp iPj H is the two-photon 



where (ev/c)Ji 
contact-interaction tensor. 

The full amplitude R = Rd + R w of inelastic Raman scattering of the photon on 
electrons changing its energy Q to Q = Q — u (u> = e/ — ti is the Raman shift), and 
momentum and polarisation from q, 1 to q, 1 is described by the Feynman diagrams 
shown in Fig. [2] The building blocks of the diagrams include Green's functions for the 
electrons and the electron-photon interaction vertices: 



G R(A) = ( £ _ H ± i0 ) 



-1 



X 
X 



evh 
Rw 



Ja-1, 



X 



evh 



(6) 
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2Vnn 

Thus, the inelastic light scattering amplitude is contributed by a one-step process R w 
(contact interaction [3H]) and a two-step process Rd involving an intermediate virtual 
state. The two-step process, such as shown in Fig. [2j consists of an absorption (emission) 
of a photon with energy Q (Q) transferring an electron with momentum p from an 
occupied state in the valence band into a virtual intermediate state, followed by another 
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electron emission (absorption) of the second photon with energy Q (fl) which moves the 
electron to the final state with momentum p + q — q, 

Rd = ^S((J q • 1)6^(^-1*) + (J-, • r)GV e/ (Jq • !))• ( 7 ) 

In graphene, the two-step process via the virtual state and the contact interaction 
process have different properties and polarisation selection rules, which reflects their 
intricate relation with the irreducible representations of symmetry group of the crystal. 
Generally, the components of the scattering amplitude realise a representation 
of the symmetry group Cq v . Since vectors 1, 1* belong to E\ representation, the R 
representation can be expanded in irreducible ones 

E 1 © E x = At © A 2 © E 2 . (8) 

Forming the corresponding combinations from the components of the polarisation vector, 
see Tabled] 

s s = |i x r| 2 , H s2 = |i-rf, s = |d| 2 = i + (ixr)(ixr), (9) 

where d = + l y I* x , lj* x - l y l*), and E s + E s2 + E Q = 2, (10) 
we can write the scattering probability in the most general form 

W = W S E S + W s2 =Ls2 + w E . (11) 

The first two terms with polarisation factor S s and S s2 , corresponding to A 2 and A\ 
representations respectively, describe the contribution of photons scattered with the 
same circular polarisation as the incoming beam. The third term, with polarisation 
factor E Q , corresponding to E 2 representation, represents the scattered photons with 
circular polarisation opposite to the incoming beam. 

Microscopically, the probability for a photon to undergo inelastic scattering from 
the state (q, q z ) with energy Q into a state (q, q z ) with energy Q, by exciting an e-h pair 
in graphene, is 

w =^sJ2 /^Pl(/l^)| 2 x/,(l-//)5( e4 + u;- e/ ), (12) 

OtiCXf 

where fi and ff are filling factors of the initial and final electronic state, respectively, 
the spin and valley degeneracies have already been taken into account, and initial 
and final states are defined by the momentum p and (sub)band index cti/f'- \i{f)) = 
|p(p + q — q, oci/f)). The angle-integrated spectral density of Raman scattering g{u) is 

•("H/ 7Sf ro ^- c ^«)' (13) 

and quantum efficiency J = J gdu expresses the total probability for single incoming 
photon to scatter inelastically on an electron and excite an electron-hole pair in the 
low-energy part of the spectrum. 
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2.1. Raman scattering in monolayer graphene 

The calculation of the Raman spectrum and intensity of the signal in each of the distinct 
polarisation components, w s , w S 2, and w Q requires the explicit expression for the Raman 
scattering amplitude R. Raman measurements are usually performed at Q, Cl ~ l-=-3 eV, 
which is much smaller than the bandwidth [17J. Raman shift, which measures the 
electronic excitation energy, is lo < 0.2 eV allowing us to expand the Green's functions 
over the large Q, and to perform summation over the intermediate virtual states of the 
process. 

Both the electron and the hole of the e-h excitation have almost the same 
momentum (p + q — q and p, respectively), since q,q C p and the momentum transfer 
from light is negligible (v/c ~ 3- 10~ 3 ), so that below we use J q = £c<5 PjP+q . For vp <^Q 
the dominant contributions to the Raman scattering amplitude are 

(14) 
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(15) 

where the main order of Rd has the matrix form of representation A 2 of Cq v . Contact 
interaction R w is responsible for the creation of the excitations with the symmetry of the 
representation E 2 , see Table [I] For scattering of photons with f2 < 70, R w <C Rd [4Tj . 

In undoped graphene the inter-band e-h pairs with ej ~ — e« ~ u/2 are the only 
allowed electronic excitations. The probability of the creating an electronic excitation 
by photons with Q < 7o, q-q<p (oj S> ^) in the main order of the expansion in 
uj/VL < 1 is 

w w he A v 2 ^- A {e s + —^—-E \ . (16) 



fi 4 1 2(6 7o r ) 

It is dominated by the contribution Rd of the first two diagrams in Fig. [2j describing 
two-step processes involving a virtual intermediate state. 

Doping in graphene, with chemical potential fi 3> flv/c blocks inter-band electronic 
excitations with energy transfer u < 2fi. After integrating over all directions of the 
propagation of scattered photons, we find the spectral density of the angle-integrated 
Raman signal, 

1 / 2 \ 2 
1„ / e v\ u) 



^~4 E '{^-c) < 17 > 
The total contribution of the electronic scattering to the total Raman efficiency is 
■n/2 / e 2 

2hc c 

In doped graphene one may also expect to see some manifestation of the intra-band 
e-h excitations in the vicinity of the Fermi level and electron energies ~ e/ ~ fi. 







g( U )du~[£--) ~10- 10 . (18) 



Their analysis requires taking into account all terms in the expression for Rd, Eqs. (15) 
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Then, we find that, for u ^ (v/c)Q <C < 70, the yield of this low-energy feature is 
61 = J 5g(u)dLu ~ KT 15 for Q ~ leV. 

Raman scattering in bilayer graphene 

For bilayer graphene with the interlayer coupling less then the base photon energy, 
71 <C f2, f2, we also study the low energy excitations in the final states, in particular, 
those with w < 71. The reason for restricting the following analysis to the low-energy 
excitations is that the higher-energy excitations decay fast due to electron-electron and 
electron-phonon processes, leaving very little chance to observe any structure in the 



Raman signal. As a result, we employ the same expansion u/Q <C 1 as in section 2.1 
Keeping terms up to the Q~ 3 order [the latter appear when the virtual state is taken to 
be in the high-energy bands] in the expansion of the Green's functions and performing 
summation over the intermediate states of the process, we get the amplitude R in the 
form of a 4 x 4 matrix 

R V " Nlxl* 

M = 





lio- x {e z x cr)a x M 
M 



Since the electronic modes with the excitation energies u < ji/2, result from the 
electron-hole pairs in the low-energy bands described by low-energy Hamiltonian 
in Eq. ([3]), we take only the part of R which acts in that two-dimensional Hilbert space, 
leaving the terms in the lowest relevant order in vp/ji <C 1 and 71 /Q <C 1. This 
projection results in effective amplitude of inelastic two-phonon process, 

2 fc2 2 

^eff ~ ^-qT- {-^(1 x 1*)* + ^ My + °M} ■ ( 19 ) 

It is important to stress that for Q ^> 71 i? e g cannot be correctly obtained within a 
effective two-bands theory, Eq. ([3]), without the described above intermediate steps, 
though for n < 71 one could formally use H q q and Eqs. ( 4|5|6|7 ) to describe inelastic 



scattering of far- infrared light. Within these approximations the scattering probability 
^ eV lL { 3< + |_=}, (20 ) 

retains the characteristic of monolayer graphene — the crossed polarisation of in/out 
photons in the dominant mode [IS]. The angle- integrated spectral density of Raman 
scattering 

^ = \{wcc)wA-- + ^ s ]^-^ < 21) 

reflects the constant density of states characteristic for the parabolic spectrum of the 
bilayers. This is different from monolayer graphene, where g(oS) oc u, reflecting the 
linear energy dependence of the density of states [47J. 
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3. Inter-Landau-level transitions in the Raman spectrum of graphene 

Upon the quantisation of electron states into Landau levels Raman scattering inevitably 
causes inter-Landau-level transitions, which, having fixed energies, manifest themselves 
on the electronic contribution to the Raman plot as a pronounced structure which can 
be used to detect their contribution experimentally [321 ED] 

3.1. Monolayer graphene 

Electronic spectrum of monolayer graphene in a strong magnetic field can be described 
as a sequence n a of Landau levels (LLs), corresponding to the states J3DJ E] 



with energy ae n , e n = y/2nhv/\B, 

where Xb = \Jh~cjeB is the magnetic length, n enumerates the Landau levels, a = + 
denotes the conduction and a = — the valence band, and $ n are the normalised 
envelopes of LL wave functions. Then, the corresponding elements of Feynman diagrams 
in Fig. [2] are 



_ qR/A _ 

e — ae n ± iO 1 
J n a n / Q ' = m#„/ jn _ie_ — c/z<5 n /__i >n e + , 

R» =^eJ-E e ±( le T)( ! * e T)- 

67 iZ ^ 

Here e± = (e x ± ie y ) is used to stress that a circularly polarised photon carries 
angular momentum m — ±1. Matrix structure of the interaction terms allows the 
following selection rules for optically active inter-LL excitations in monolayer graphene: 

i) n-^n + ; ... n~ -> (n + 1) + , , . 

U) (nTl)^(n±l) + ; mJ (n + 1)" -> n+. 1 J 

The excitation of the e-h pairs by Raman scattering in graphene at strong magnetic fields 
produces the electronic transition i) between LLs, with angular momentum transfer 
Am = and excitation energy u = 2e n , see Fig. [3j and transitions ii), with Am = ±2 
and oj = e n -i + e n +i. The amplitudes of these two processes, 

l(e^) 2 ^[ (le+)(l*e_) (le+)(l*e_) 



a=± 



ft — e n — ae n+ ± e n - Q - ae n ^i 
(le_)(l*e + ) + (le_)(l*e + ) 



— e n — ae n+ \ e n — Vl — ae n -\ 

1 (evhf 



R(nTi)-+(n±i)+ = =F j-^-(le±)(l*e±) x 
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Table 2. Raman-active inter-LL excitations in graphene. 
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^ — £n+i — QJ£ n £ n -i — ^ — Q!£ n 



are such that R n -^ n + 3> -R( n =Fi)-->(n±i)+ f° r oj <^ VL. The latter relation is determined 
by a partial cancellation of the two diagrams constituting R^. Notice that the inter-LL 
modes n~ — > n + have the symmetry of the representation A 2 in Table [2] and the same 
circular polarisation of "in" and "out" photons involved in its excitation. Finally, the 
contact term R w in Fig. [2] allows for a weak transition iii) n~ — > (n ± 1) + , with the 
amplitude R w <C R n -^ n +. This transition, together with selection rules Am = ±1 
resembles the inter-LL transition involved in the far-infrared (FIR) absorption [21 E]. 
However, the FIR-active excitation is 'valley-symmetric' [51] and belongs to the 
representation E 1; whereas the Raman-active n~ — > (n ± 1) + mode belongs to E 2 , 
allowing the latter to couple to the T-point optical phonon and, thus, leading to the 
magneto-phonon resonance feature in the Raman spectrum [29]. Also, R w originates 
from the trigonal warping term in % which violates the rotational symmetry of the 
Dirac Hamiltonian by transferring angular momentum ±3 from electrons to the lattice, 
allowing thus the change of angular momentum by ±1. 

The dominant inter-LL transitions n~ — > n + determine the spectral density of light 
scattered from electronic excitations in graphene at high magnetic fields: 

(2 2 / \ \ 2 
^ ( ~^f~) ~2e n ). (24) 

' n>l 

Here we use Lorenzian j n (x) = •n~ 1 Y n /[x 2 + T 2 ], and inelastic LL broadening T n which 
increases with the LL number, and the factor S s in Eq. (24) indicates that "in" and 
"out" photons have the same circular polarisation. 



3.2. Bilayer graphene 

Same as in non-magnetic case, e-h excitations in bilayer graphene are created in low- 
energy Landau levels, as at high energies the Landau level broadening due to, for 
example, electron-phonon interaction, will smear out the LL spectrum. In strong 
magnetic fields, low-energy Landau levels in bilayer graphene can be described [52] 
by 



|n ° ) = ^UI-J forn - 2 ' |0/1> = l >• (25) 
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a) 



b) 




i 2 + 



] 1 + 





Figure 3. Spectral density g(uj) of light inelastically scattered from electronic 
excitations in undoped monolayer (a) and bilayer (b) graphene at quantising magnetic 
helds (solid line) and in doped graphene at B — (dashed line). Sketch illustrates 
selection rules for processes corresponding to the visible features at Raman plot. e n 
denotes the LL energies in graphene. In addition to dominant n~ — > n + processes 
at both plots, in bilayer spectrum the feature corresponding to 2~ — > and — !• 2 + 
transitions can be recognised. 



h 2 v 2 , 

with energy e n a = 2a — r^-yn(n — 1). 

In a neutral bilayer, all LLs have additional four-fold degeneracy (two due to the electron 
spin and two due to the valley). The peculiar n = and n — 1 LLs are degenerate 
at e = producing to an 8-fold degenerate level. Projecting the effective transition 
amplitude R e Q onto the eigenstates \n a ) we find the electronic Raman spectrum. This 



procedure leads to the same selection rules as in monolayer graphene, see Eqs. (23). 
Among those, i) are the dominant transitions. For a undoped bilayer, the angle- 
integrated spectral density g(ou) of Raman scattering in the strong magnetic field is 
equal to: 



v 7 k n>2 



nJ 



X 
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+ 7,^2^-^+1)+ - e ("-l)+)) f- ( 26 ) 



* n>3 

The first term describes the dominant contribution due to the n~ — > n + transitions, 
while the last two terms describe the spectral density of the (n ^ 1)~ — > (n±l) + 
transitions. 

Low-energy electronic contribution to the Raman plot is shown with a solid line 
in Fig. [3]d. The dominant peaks are due to the n~ — > n + transitions. The quantum 
efficiency of a single n~ — > n + peak in Fig. |3Jd is similar to value I\ in monolayer 
graphene [17] . 

A weaker feature in Fig. [3|d corresponds to both 2~ — > and — > 2 + transitions 
and is visible as it is positioned to the left of the first 2~ — > 2 + peak. The quantum 
efficiencies of the (n±l)~ — > (n=)=l) + transitions are smaller by the factor (32~) m 
comparison to the n~ — > n + transitions. This is different from the monolayer graphene 
case, where the corresponding ratio between quantum efficiencies of (n±l)~ — > (n=pl) + 
and n~ — > n + transitions is (^) , much smaller than for the bilayer. 

4. Magneto-phonon resonance in homogeneously and inhomogeneously 
doped and strained graphene 

Single- and multiple-phonon-emission lines in the Raman spectrum of graphene and 
the influence of the electron-phonon coupling on the phonon spectrum have been 
investigated in great detail [211 [25l [26l EZl [531 EH ESI EOl EU [32]. The coupling of 
electromagnetic field and phonon excitations is realised via the excitation of the virtual 
e-h pair, with the same symmetry properties as the T-point optical phonon which is 
created upon the recombination of the pair. Analysing all possible combinations of 
polarisation vectors one may see from Table [T] that the amplitude of Raman scattering 
with excitation of the phonon mode is proportional to 



that describes quantised sublattice displacement 



U(r, t) = (U A - U B )/V2 = T J^7, { u ^e l(kr ^ ot)/h + h.c.) , (2* 



where M is the mass of single carbon atom, u is a polarisation of the phonon, and 
expression for Rq can be taken from Basko (53]. Raman efficiency then is proportional 
to 

tfoc^dXImLd, (29) 

vv 

where D is the phonon propagator renormalized by electron-phonon interaction. The 
latter strongly affects the energy of the phonon under the conditions of its resonance 
mixing with the inter-LL electronic transitions, the effect known as magneto-phonon 
resonance. 
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We use both phenomenological approach and the tight-binding model to study the 
interaction of electrons with T-point optical phonons in case of strong magnetic field. 
The phenomenology is suitable to describe the low-energy excitations of electrons in 
the vicinity of Dirac points (K-points), while tight-binding model is needed to perform 
the integration over the whole Brillouin zone [53] when estimating the intensity of the 
G-phonon-line in Raman. To discuss the electron-phonon interaction it is instructive to 
trace the electron Hamiltonian 

H = H - Q • U (30) 

back to the tight-binding model 

where k is a momentum within the Brillouin zone, K = (— 0), v = ^joa/h, vectors 
ri = (— f , ^§)j r 2 = (f , ^tj)j an d 1*3 = (0, — -j=) are the positions of nearest "B" atoms 
with respect to "A" atom, and F = ^= where ^ is a response of the hopping 
amplitude to the variation of the bond length. 

The renormalisation of the phonon energy by its interaction with electrons is taken 
into account in the random phase approximation. Defining electrons' bare Green's 
functions phonon bare propagator, and electron-phonon interaction vertex, 

G R/A = (e-H±iO)-\ 
G K = {l-2f){G R -G A ), 
_ 2u 




Mu 

we use Dyson equation to find the renormalised phonon propagator: 

x Tr{Q v G R (e)Q u G K (e + u) + Q V G K \e)Q v G A {e + u)}, (33) 
D = === = — + (03 + cQd cQd +... (34) 



2 



— o== 



Here A = (f^) 2 = ^ (jt^dr) ~ 0,018 is dimensionless constant of the electron- 
phonon coupling (26]. Trace is taken in AB sites space and a summation over spins is 
already performed. As momentum A; of a phonon created by the Raman process is small 
by parameter v/c, in the following it can be neglected. 
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The integral for polarisation loop in Eq. (33) is divergent in the frame of linear 
Dirac model, so it can be split into two parts: 



Un5 w + An, 



n 



-0.64 7( A 



(35) 



Here Ho is calculated numerically within tight-binding model, and /SXl vv is determined 
within the effective Dirac model, after subtracting the divergent part from the integrand. 
For B = the answer is An^, oc 5 VV and it repeats the derivation of Raman line shape 
found in Ref. [281 [2S1 EU H |E2 HE [32] 



9 



g: 



(w-w g ) 2 + 7g' 

where the width of the G-peak and its position are determined by 

A 



7g = l P h ~ 7W0 - 2/i), 



00 G = w 



4 

A^ + A^ 

7T 4:71 



log 



2/i + oo 



(36) 
(37) 



2/i — uj 

where 6o = oo + 0.64A7 (we implied that oo ^> 0.64A7 ^> Xep), and expression for I G 
can be taken from Basko 1531. 



4-1. Magneto-phonon resonance in a homogeneously doped monolayer 

In the presence of magnetic field the continuous electronic spectra splits into the set of 
discrete Landau levels (LL). Working in the basis of LLs Eqs. ( [22|33|34] ) and choosing 
the circularly polarized basis e / = ^(e^ ± ie y ) for incoming and scattered light, we 
find that 



(e )„(e )*n + (e )„(e )*n , where II 



o/o 



n + An 0/0 . 



For the chemical potential pinned at the zero LL 



An 



o/o 



c-2 

zl 

71 



-fo l-/o 1 
—— + — — + u 

Ei±U —Si ±00 E\ 



E 

n>2 



(38) 



where fo is a density-dependent filling factor of the zero level and fl n ,± = £ n ± £n-i 
are the frequencies of electron- hole excitations. If /i is large enough, so that m level 
(m > 1) is partially filled, then 



An 



0/0 



e-2 

2tt 



fr, 



-fr, 



+ 



fr, 



Q m _ ± 00 Vl m +i^- ± 00 Q-m + ± 00 



l-fn 



1 



-n 



m+l,+ 



- 2oo 2 

4 



E 



1 



.(39) 



n>m+2 

Then, the net Raman efficiency of the T-point phonon mode in graphene in a strong 
magnetic field is 



9 



E i de 

.4=0/0 



[00 



«„-AAn»)H 7p 2 ; 



(40) 




Figure 4. G-peak splitting in magnetic field. Density plots represent the dependence 
of the Raman efficiency on Raman shift and magnetic field, a) Schematic plot shows 
phonon line, dashed curves represents dominant Raman active transitions nT — > n + , 
which do not interact with phonon line, and solid curves are weakly coupled Raman 
active transition nT — > (n±l) + , which interact with phonons forming the anticrossings; 
b) first anticrossing due to the transitions 1~ — » 0, — > 1 + ; c) 2 nd , 3 rd , and part of 
4 th anticrossings. In b) and c) plots are drawn for a zero electron density, so that the 
contributions from different polarisations are equal. 



where the two terms in the sum correspond to the opposite circular polarisations of the 
optical phonon, and it is plotted in Fig. [4] for several doping densities. 

The line described by the Eq. (40 ) displays two types of typical behaviour depending 
on the range of parameters. Far from the resonance only a slight shift of the electron 
and phonon excitations exists, 

w 0/0 Riw + AAn 0/0 (uo). (41) 

Since the line shift may differ for the opposite circular polarisations of the phonon, this 
leads to a small splitting of the G-line in Raman. In the vicinity of the intersection of 
the phonon line with the electron-hole excitation line for several doping densities the 
polarisation operator becomes large and causes significant change in spectra, forming 
an anticrossing, in addition to the line splitting. The form of the line is determined by 
four solutions 0Jq/ o 

w-wo-AAn 0/0 ( W ) = 0, (42) 
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two for each polarisation with energies above (+) or below (— ) the bare phonon line. 
The form of such anticrossings in undoped graphene is shown at Fig. |4| The resonance 
takes more sophisticated form in the case of the doped graphene. With increasing of 
the electron density (and, correspondingly, the filling factor), the transition 1~ — > + 
invoked by anticlockwise polarised light is weaken by Pauli blocking (see first anticrossing 
at Fig. [5^,, while the CT — > 1 + transition remains. When the n = Landau level 
is fully occupied (n e = 1/tt\ 2 b w 1.4 x 10 12 cm -12 ,), in anticlockwise polarisation the 
anticrossing completely blocked (see Fig. [5]d, which coincides with the experimentally 
observed features [4"9] . 

4-2. Influence of inhomogeneous doping and strain on the magneto-phonon resonance 
lines shape 

In graphene deposited on the substrate, the density of the electrons is very often, 
not the same everywhere in the sample, but varies, forming puddles at the length of 
100 -T- 500nm (56]. Also, ripples of graphene sheet, induce inhomogeneous strain. At 
weak magnetic fields the the electron density and strain disorder lead to the scattering 
of the charge carriers. The situation changes when the external magnetic field (~ 25T) 
is applied. At such a strong field, the magnetic length is about 5nm, which is much 
smaller than the typical size of the spatial fluctuations of both carriers density and 
strain. In this case, the Raman response of the graphene flake is the average of the 
locally defined spectra. Below, we discuss how such averaging changes the fine structure 
of the magneto-phonon resonance. 

Inhomogeneous density The spatial distribution of the electrons density, results in 
summing up of the local spectra contributions from the various parts of the graphene 
flake. The net picture is the overlay of the anticrossings at different filling factors (like 
one shown at Fig. |5^i,b) weighted with the normal distribution factor, which we use 
mostly for the illustration purposes. 

As shown in Fig. [5}:,d, this leads to the inhomogeneous broadening of the side 
lines related to the variation of the size of the anticrossing gaps between the - — » 1 + 
transition and the phonon, and it blurs out the much weaker splitting/anticrossing of 
the phonon and 1~ — > + transition. 

Inhomogeneously strained graphene Inhomogeneous strain in graphene is often 
associated with magnetic field, -Entrain for electrons in K valley and — -Bstram in K', 
which estimated as -B s tram ~ where a is the bond length, I and h are the linear 

scale and height of a ripple in graphene [31] . In the external magnetic field B the total 
effective field acting on electrons in the valleys K and K' is different, 

B^ = B + £-B s train, 

where £ = + for K point and £ = — for K' . Difference in magnetic field leads to the 
different Landau spectra in each valley, thus, splitting the values of the external field B, 
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Figure 5. Anticrossing of 1~ — > + and 0~ — > 1 + transition with the optical 
phonon for various levels of graphene doping, a) homogeneous graphene, electron 
density n e = 1 X 10 12 cm~ 2 : transitions in both circular polarisations are active, 
but l - — > + is suppressed; b) graphene with a homogeneous doping with electron 
density n e — 2 x 10 12 cm -2 : only CP — s- 1 + transition is active; c) graphene with 
inhomogeneous electron density, such that the average density is n e = 1.5 x 10 12 cm~ 2 
and distribution with the variance Sn e — 0.4 x 10 12 cm~ 2 . d) The same as in c), but 
with n„ = 2.3 x 10 12 cm -2 and 8n e = 0.4 x 10 12 cm -2 . 
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Figure 6. The CP — > 1 + and 1~ — > + anticrossing for strained graphene with 
electron density n e = 2 x 10 12 cm~ 2 . a) homogeneously strained graphene with 
effective -Bstrain = 3T: one anticrossing splits into two shifted vertically by 2i? stra in; b) 
graphene subjected to inhomogeneous strain averaged over the distribution with the 
variance 5i? s train = 3 T: the peak corresponding to 0~ — > 1 + transition has a distinctive 
asymmetric shape. 



for which, e.g. the resonance of K — >■ 1^ and K , — > \\, with the phonon is realised. 
Formally, this can be taken into account by a substitution AU(B) — > | ^5=± 
in Eq. Q. 

For the local Raman spectrum related to an area with a fixed B strain , this results 
in a fine structure of the anticrossing shown in Fig. [6^,, where an additional line appears 
seemingly crossing the phonon mode without splitting. Such a behaviour occurs due 
to the fact that instead of an anticrossing/mixing of two modes, now, three modes are 
coupled, forming one phonon-electron mixed excitation with energy pinned in the middle 
between the two strongly split side bands. Here, the central phonon line corresponds 
to the circular polarisation of the phonon for which mixing with electronic modes is 
suppressed by Pauli blocking, due to high density of electrons used in the calculation. 
Averaged out a normal distribution of -Bstrain with zero average expectation and variance 
<5-£>strain = 3T, the anticrossing branches become asymmetrically broadened, as shown in 
Fig. §3. 



5. Summary 

We presented a theory of inelastic scattering of photons in mono- and bilayer graphene 
accompanied by the excitation of electron-hole pairs. The dominant scattering processes 
are characterised by the crossed polarisation of in/out photons while the Raman 
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efficiency is proportional to the density of states (linear increase with Raman shift for 
monolayer and constant for bilayer) and in case of doping the gap of width 2[i opens 
up. 

In strong magnetic field Raman spectra acquires the pronounced structure with 
peaks corresponding to the inter-Landau-level transitions. The selection rules for the 
strongest mode preserve LL number n: n~ — > n + , and quantum efficiency of the single 



peak in Fig. 



is Ji ~ e Isf ) ~ 1CT 12 per incoming photon for Q ~ leV photons 
in magnetic field B ~ 10T. The theoretical prediction for the intensity of the phonon- 
induced G peak, [53] which is a well known Raman feature in carbon materials |54j, 
estimates Iq ~ 1CT 11 . This result is close in the order of magnitude to the intensity of 
a single n~ — > n + peak. So, the spectral features of inter-Landau-level transitions in 
graphene can be observed experimentally [H]. The weaker mode with same selection 
rules, similar to the absorption process, n~ — > (n ± 1) + , but belonging to different 
symmetry group representation (E2, while absorption has Ei), is coupled to the phonon 
mode, creating the anticrossing and the G-peak line. 
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